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ABSTRACT 

Accretion from disks onto young stars is thought to follow magnetic field 
lines from the inner disk edge to the stellar surface. The accretion flow thus 
depends on the geometry of the magnetic field. This paper extends previous 
work by constructing a collection of orthogonal coordinate systems, including the 
corresponding differential operators, where one coordinate traces the magnetic 
field lines. This formalism allows for an (essentially) analytic description of the 
geometry and the conditions required for the flow to pass through sonic points. 
Using this approach, we revisit the problem of magnetically controlled accretion 
flow in a dipole geometry, and then generalize the treatment to consider magnetic 
fields with multiple components, including dipole, octupole, and split monopole 
contributions. This approach can be generalized further to consider more complex 
magnetic field configurations. Observations indicate that accreting young stars 
have substantial dipole and octupole components, and that accretion flow is 
transonic. If the effective equation of state for the fluid is too stiff, however, 
the flow cannot pass smoothly through the sonic points in steady state. For a 
multipole field of order £, we derive a general constraint on the polytropic index, 
n > i + 3/2, required for steady transonic flow to reach free-fall velocities. For 
octupole fields, inferred on surfaces of T Tauri stars, the index n > 9/2, so that 
the flow must be close to isothermal. The inclusion of octupole field components 
produces higher densities at the stellar surface and smaller areas for the hot spots, 
which occur at higher latitudes; the magnetic truncation radius is smaller (larger) 
for octupole components that are aligned (anti-aligned) with the stellar dipole. 
This contribution thus increases our understanding of magnetically controlled 
accretion for young stellar objects and can be applied to a variety of additional 
astrophysical problems. 
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Introduction 



During the star formation process, most of the mass that becomes part of the nascent 
star initially falls onto a surrounding circumstellar disk, rather than directly onto the stellar 
surface. This accretion process can be (indirectly) observed during the T Tauri phase of 
evolution, after the star/disk system has emerged from its protostellar envelope. The transfer 
of material from the disk to the star takes place at the inner disk edge, which generally does 
not extend to the stellar surface. Instead, the inner edge is connected to the star through 
magnetic fields, and accretion takes place along these field lines. 

This paradi gm of rnagne t ic accretion was developed in the earl y 1990s for classical T 



Tauri stars (e.g., Konigl 1991 ; Shu et al. 1994 : Hartmann et al. 1993) and is roughly analo- 



gous to that of magnetically controlled accretion from disks onto neutron stars (iGhosh fc Lamb 
19781 ) and black holes (Blanford & Payne 1982; see Uzdensky 2005 for further discussion). 
Many of the diverse observational characteristics of accreting T Tauri stars can be explained 
within the basic framework of the magnetospheric accretion scenario. The shapes of spectral 
energy dist ributions (SEDs) in the near infrared are cons i stent with cavities in inner dust 
disks (e.g., iKenyon fc Hartmannlll987l: lAdams et al.lll988l : iRobitaille et al.l 120071 ). although 
gas likely extends closer to the star ( iNajita et al.ll2003l ). The excess emission, primarily at 
IR and UV wavelengths and apparent from the SEDs of accreting T Tauri stars, can be 
explained by the reprocessing of the stellar photons by dusty material in circumstellar disks 
and from shock emission at the base of the accretion columns. Gas in the accretion flow 
rains down onto the stellar surface p roducing hotspots that radiate primarily in the UV 



and also in the soft X-ray waveband (IKastner et al. 



20021 : lArriroffi et al.l 1201 If ). Accretion 



related hotspots, in addition to cool spots which arise where bundles of magnetic flux rise 
through the stellar surface in to the atmosphere, contribute to the high level of photometric 
variability of T Tauri stars (IBouvier et al.lll993[ ). Spectroscopically, the photospheric ab- 
sorption lines of accreting T Tauri stars are typically shallower than those observed from 
non- accreting T Tauri sta rs and main-sequence stars of the same spectral type (see Figure 
1 of iHartigan et al.lll991l). as a result of the additional continuum emission from the ac- 
cretion spots (jCalvet fc Gullbringlll998l ). Many emission lines often exhibit red-shifted and 
blue-shifted absorption con iponents, sometimes simultaneously, cha racteristic of accretion 
and outflows, respectively (lEdwards et al.l Il994l : iFischer et al.ll2008l ). and suggesting that 
the star-disk interaction region contains complex kinematic gas flows with material both 
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accre ting; onto the star and being launched from the system in outflows (see iBouvier et al. 
2007 for a review). 



A key assumption of the magnetospheric accretion model is that T Tauri stars sup- 
port large-scale magnetic fields that are sufficiently globally ordered and strong e nough to 
disrupt the disk at a distance of few stellar radii. As demonstrated by iKonigll ( 1199 ll ). a 
stellar dipole field of polar strength B ~ 10^ G is sufficient. Init i al magnetospheric accre- 
tion models focused on dipole magnetic fields (see also Lil 19961: Li fc Wilson 1999), and 
have included detailed heating an d cooling calculations (Martin 1996h. polytropic equa 



Zanni &: Ferreira 



axis (IRomanova et al 



tions of state (|Koldoba et al.ll2002l ). detailed numerical treatments (IRomanova et al.l 12002 



20091). and dipole fields that are tilted with respect to the stellar rotation 



2003|). 



Strong stellar-disk-averaged surface fields have now been measured on a number of T 
Tauri stars in different star forming regions, most successfully through the detailed anal- 
ysis of magnetically sensitive, and therefore Z eeman broadened, lines in intensity spectra 
(jJohns-KrullI 120071 : lYang fc Johns-Krulll 1201 ll ). However, such broadening measurements 
give no information about the stellar magnetic field topology. Another manifestation of the 
Zeeman effect, namely the circular polarization of magnetically sensitive lines, does yield in- 
formation about the field geometry. Initial spectropolarimetric studies show mixed results: 
Attempts at measuring the polarization signal in ph otospheric absorption lines often fail to 



detect the presence of surface magnetic fields (e.g., iJohnstone fc PenstonI Il986l ). However, 
as opposite polarity (positive and negative) surface field regions give rise to signals that are 
polarized in the opposite sense, a net circular polarization signal of zero i s consistent with 
T Tauri stars hosting complex surface fields (IValenti fc Johns-Krulll |2004| ). In contrast, a 
strong and rotationally modulated circular polar ization signal was detected in the Hel D3 
(5876A) emission line by I Johns-KruU et al.l ( Il999l ). This particular line of helium has a high 
excitation potential and is thought to form at the base of accretion columns. The strong 
rotationally modulated signal in this, and other accretion related emission lines, is found to 
be well described by a simple model where t he bulk of the accreting gas lands on the stellar 
surface in a single polarity radial field spot ( IValenti fc Johns-Krulll l2004l ). This finding sug- 
gests that even though accreting T Tauri stars host complex surface magnetic fields, their 
large-scale field topology, and in particular the portion of the field that carries gas from the 
inner disk to the star, is simpler and globally well-ordered. 

The non-dipolar nature of T Tauri magnetic fields has recently been confirmed. Spec- 
tropolarimetric observations, combined with tomographic imaging techniques whereby the 
rotational modulation of the Zeeman signal is model ed, have allowed magnetic surface maps 



to be derived for a number of accreting T Tauri stars ( iDonati et al.ll2007l l2008l : iHussain et al. 
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20091 : iDonati et al.ll2010al Jbl. 12011 



For completeness we note th at magnetic map s have 



also been pubhshed for a handful of non-accreting T T auri stars (IDunstone et al.l 12008 



Skelly et al.l l2010l : iMarsden et al.l 1201 ll : IWaite et al.l l201ll ) . These maps are constructed by 
considering the rotational modulation of the polarization signal detected in both photospheric 
absorption lines, which form uniformly across the entire stellar surface, and the signal in ac- 
cretion related emission lines, which trace the field at the base of accretion columns. In 
practice, the polarization signals detected in photo spheric absorption lines are small, and 
cross-correlation techniques (e.g.. iDonati et al.lll997l) are employed in order to e xtract infor- 
mation from as many spectral lines as possible (see iDonati &: Landstreetl |2009| for a review 
of Zeeman-Doppler imaging, the technique used to construct stellar magnetic maps, and 
Donati et al.l (|2010bl ) for its specific application to T Tauri stars). 



The observationally derived magnetic maps can be decomposed into the various spherical 
harmonic modes. Some stars are found to host very complex magnetic fields with many 
high order components and strong toroidal field componer its, such as both st ars of the 
close binary V 4046 Ser dPonati et all boilch. V2247 Oph JPonati et al.lboiOah . CR Cha 



and CV Cha (iHussain et al. 



200 



Intriguingly, however, many accreting T Tauri stars 
have field topologies that are well described as dipole-octupole composite fields. The polar 
strength of the dipole and octupole field component varies from star t o star. AA Tau has a 
dominantly dipolar magnetic field, with a weak octupole component ( IDonati et al.ll2010bl ). 
whereas TW Hya hos ts a dominantly octupolar magnetic field with a weak dipole component 
( IDonati et al.ll2011bl ). The same is true for V2129 Oph, although the dipole component has 
been observed to vary by a factor of three (from a polar strength of ~ 0.3 to ~ 0.9 kG) over 



a time scale of four years perhaps hinting at the existence of a magnetic cycle (IDonati et al. 



2011al ). BP Tau, one of the best studied accreting T Tauri stars, hosts a magnetic field with 
both strong dipole and octupole field components (IDonati et al.ll2008l ). 



These observations provide clear motivation to consider magnetic fields with both dipole 
and octupole components in models of accretion flow, and motivated (at least in part) by 
the availability of the new observational data, magnetosphe ric accretion raodels with higher 



order multi p ole stellar magnetic fie l ds have been develop ed ( IGregory et al.ll2006t iLong et al. 
2OO7I . boosi : iMohantv fc SM boosi : ICregorv et aPboiol ). To date, flow models with ac- 



cretion taking place along dipole stellar field lines have successfully reproduced emission 



line profiles and their rotational yariabi lity, including helium lines (iBeristain et al.l 12001 



Fischer et al. 



2008 



lines (Hartmann et a. 



1994 



Kurosawa et al.ll2011l). calciurn lines (lAzevedo et al.ll2006l) and hydrogen 



MuzeroUe et al.ll2001t ISymington et al. 



2005 



Lima et al. 2010 



Kurosawa et al.ll201ll ). However, the simulated line profiles based on dipo le stellar mag- 
netospheres show more variability than is observed (jSymington et al.ll2005[ ). and for some 
stars the inferred trajectory of accretion flows close to the stellar surface, where the higher 
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order field components will influence the infalling columns o f gas ( iGregory et al.ll2008l ). are 
inconsistent with the dipole flow model ( iFischer et al.ll2008l ). 



Deep, broad, and often rotationally modulated red-s hifted absorption coniponents are 
commonly detected in accretion related emissio n lines (e.g. jEdwards et al. I ll994i iBouvier et al, 
20031 : lEdwards et al.ll2006l : iFischer et al.ll2008l ). Analysis of such features allows the determi- 
nation of the kinematic gas flow (the accretion column) crossing the line-of-sight to the star. 
The material arrives at the stellar surface with highly supersonic speeds of several hundred 
km s~^. Since the inflowing material leaves the inner disk edge at low (subsonic) speeds, the 
gas must make a sonic transition while it follows magnetic field lines from the disk onto the 
star. As a result, accretion is described by transonic flow soluti ons (which a re essentially 
the reverse of the well-known Parker model of the Solar wind; see lParkerlll965[ ). As we show 
in this paper, the requirement of a smooth transition through the sonic point, the need for 
free-fall speeds in the inner limit, and the divergence properties of higher order multipole 
moments, jointly place strong constraints on the allowed polytropic index of the accretion 
flow. 

The goal of this paper is relatively modest: Building on the theoretical work outlined 
above, this paper provides an analytic, or at least semi- analytic, treatment of magneti- 
cally controlled accretion flows (where the term "semi-analytic" refers to models where the 
equations are reduced to, at most, ordinary differential equations). Motivated by the obser- 
vational work outlined above, we focus on transonic solutions and generalize existing work 
to include higher order multipoles, especially magnetic fields with both dipole and octupole 
components. To reach this objective, we construct (novel) orthogonal coordinate systems for 
each given magnetic field configuration under consideration. One of the coordinates (denoted 
here as p) follows the magnetic field lines, whereas the other coordinate (denoted here as q) 
is orthogonal to the first in the poloidal plane (see Figured]). We specialize to the case of 
axial symmetry, which applies to systems where the magnetic field configuration co-rotates 
with the central star, and where the field is strong (so that the toroidal field component is 
small and the angular velocity along a streamline is constant). 

Because this paper constructs coordinate systems {p, q) in the poloidal plane, it is im- 
portant to outline why this framework is useful: The coordinate p follows the magnetic field 
lines, and hence the streamlines, so that the fluid fields are functions of the coordinate p, and 
the value of p measures the position along the field line. By construction, Vp is parallel to the 
magnetic field. As we show below, the divergence operator that describes the flow includes 
the quantity |Vp|, which is the inverse of one scale factor of the coordinate system (e.g., 
Weinreich 1998). This scale factor must be included to properly describe the divergence and 
hence the physics of the flow. The perpendicular coordinate q labels the streamlines, i.e., the 
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flow follows lines of constant q. The flow trajectory is thus specified by the functions ^(5) or 
equivalently 9{^) along a streamline and 9 are spherical coordinates); these relations are 
necessary to evaluate the relevant functions along the field lines. As a result, one must have 
a description of the scalar fields p and q that make up the coordinate system in order to solve 
for the flow. However, one could, in principle, include the correct scale factor {hp = |Vp|~^) 
and use the functions ^(6*) that specify streamhnes (lines of constant q) without an explicit 
reference to the fact that p and q are coordinates. This strategy has been used previously for 
dipole fields (Blandford & Payne 1982, Hartmann et al. 1994, and others). For more com- 
plicated magnetic field configurations, considered here, it is more straightforward to define 
the coordinate system and work within it — we gain additional physical understanding by 
being aware of the coordinate system in which the flow is taking place (for dipole fields, such 
coordinates have been used to study accretion onto white dwarfs - see Canalle et al. 2005, 
Saxton et al. 2007). Finally, the resulting framework can be readily generalized for more 
complex magnetic fields (Appendix |A]) and can be used in a variety of other applications 
(e.g., outflows from Hot Jupiters, Adams 2011). 

This paper is organized as follows. We start with a general discussion of coordinate 
systems in Section [21 The equations of motion for fluid flow are then considered in Section 
[3l where they are formulated in terms of these new coordinate systems; by construction, the 
fluid fields are functions of only one variable, which measures the location along the magnetic 
field line(s). We then derive a general constraint on the polytropic index n required for the 
flow to reach free-fall speeds; since n must be large, the flow can be described to good 
approximation using an isothermal equation of state. In Section |U we consider the magnetic 
field to be that of a pure dipole; although this case has been studied previously, we re- 
formulate the problem in new coordinates to illustrate our approach. In Section [5], we 
consider more complicated magnetic field configurations including both dipole and octupole 
components (consistent with current observations of accreting young stars). The inner disk 
edge is truncated magnetically, and hence the inner boundary depends on the magnetic field 
structure. Using the magnetic field configurations of this paper, we revisit the magnetic 
truncation radius in Section [61 The paper concludes, in Section [71 with a summary and 
discussion of our results. During the course of this work, we have derived a number of 
mathematical results that apply to general coordinate systems of the form considered herein; 
these results are collected in Appendix [XI For completeness, we also consider magnetic fields 
with dipole and radial (split-monople) components (Appendix [B|); this configuration arises 
when the inf all- collapse flow that forms the disk drags in magnetic field lines from the original 
molecular cloud core, and also when stellar winds open up field lines to become (nearly) 
radial. Finally, Appendix [O provides a consistency check by showing how this formalism 
explicitly conserves mass. 



X 



Fig. 1. — Schematic diagram showing the coordinates used here, for one quadrant of the 
poloidal plane. The star is shown in the lower left corner. The curve depicts one magnetic 
field line that connects the stellar surface to the inner edge of the disk. At each point along 
the field line, basis vectors describe the coordinate system. At one such point, the basis 
vectors are shown. The field fines correspond to fines of constant q, so the basis vector = 
Vq is perpendicular to the field line. The coordinate p provides a measure of the distance 
along the field line, so that the basis vector = Vp points along the field line. 
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2. Construction of Coordinate Systems 



In this treatment, we develop a series of coordinate systems tliat follow the magnetic 
field lines. We thus start with a given magnetic field configuration B{r,6), which we take 
to be axisymmetric (the magnetic field is poloidal, with no toroidal component, = 0). 
The origin is located at the center of the star. Since the star is expected to be rotating, 
the coordinate system co-rotates with the star, and is thus fixed. The effects of rotation 
are then incorporated by including (the usual) non-inertial terms in the equations of motion 
(see Section [3]). By definition, coordinate systems are defined by scalar fields, which we 
denote here as {p, q). If the magnetic field is curl- free (and hence current-free), then the field 
B can be written as the gradient of a scalar field. The coordinate p specifies the distance 
along the field line and hence the gradient Vp points in the direction of the field line. The 
coordinate q is perpendicular so that Vg ■ Vp = 0. Notice that lines of constant q correspond 
to the magnetic field lines; similarly, lines of constant p correspond to equipotentials of the 
multipole field under consideration. The coordinates (p, q) are thus orthogonal coordinates 
in the poloidal plane, and can be used instead of the more familiar spherical coordinates 
(r, ^). The azimuthal angle is the same for both cases and provides the third orthogonal 
coordinate. We could also use cartesian coordinates, where the pair (x, z) define the poloidal 
plane. 

As outlined above, this treatment neglects the azimuthal field component. In some 
systems, the presence of a twisted field component {B^ ^ 0) above the disk is necessary 
for maintaining balance in the angular momentum transport (Long et al. 2005). However, 
recent numerical simulations (e.g., Romanova et al. 2011, Long et al. 2011) indicate that 
the initial potential field structure is retained for regions within the disk truncation radius, 
so that our approximation [B^ = 0) is viable. For completeness we note that, in general, 
the field will not strictly be a potential field when the azimuthal component is nonzero 
(although a potential field provides a good approximation); the azimuthal component also 
affects angular momentum transport. 

Specification of a coordinate system requires not only the coordinates themselves, but 
also the basis vectors. For orthogonal coordinate systems, one can also use the scale factors. 
The set of covariant basis vectors e^- arises from the gradients of the scalar fields that define 
the coordinates (Weinreich 1998), in this case = Vp and = Vg. Note that these 
quantities are basis vectors, rather than unit vectors, so that their length is not, in general, 
equal to unity. The corresponding scale factors hj are given by the relation 



The schematic diagram of Figure [T] depicts one magnetic field line and the basis vectors at 
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one point along the field line. The coordinate p measures the distance along the field line 
and the basis vector = Vp points in the direction of the field line. The coordinate q is 
constant along field lines, and thus labels both the field lines and the streamlines. The basis 
vector = Vg is perpendicular to the field line. 

In magnetically controlled accretion, one important geometrical effect is that the diver- 
gence operator takes a non-standard form. For orthogonal coordinate systems, the general 
form of the divergence operator is given by 



V- V 



hphqh^ 



d d 

^ (hgh^Vp) + — (hph^Vg) 



+ ^^. (2) 
^sin^ 90 ^ ' 



The quantities hp, hg, and are the scale factors of the coordinate system, as defined by 
equation ([T]). 

In this application, we consider the fields to be axisymmetric so that the derivatives 
vanish. Further, for fiow along field lines, the vector fields (e.g., the velocity field) have only 
one component and depend on only one coordinate, so that the divergence operator collapses 
to the form Id 1 dV V d 

^ hphgh^ dp^ hp dp ^ hphgh^ dp^ '^^ ^ ^ 

As shown in Appendix lAl hqh^ = hp and hp oc |B|~^ (see Result 1 and the text below). 
For flow along the field lines, the divergence operator thus has a form consistent with the 
well-known result that the effective area element of a moving fluid element varies inversely 
with the fleld strength (Parks 2004). 



3. Fluid Equations for Magnetically Controlled Flow 

In this treatment, we consider steady-state solutions and assume that the magnetic fleld 
structure is strong enough to dominate the flow. As a result, the magnetic fleld is held flxed 
and fluid elements must follow the fleld lines. The magnetic fleld is assumed to rotate with 
the star, so we work in a rotating reference frame with rotation rate fl = Qz. Under these 
approximations, the equations of continuity, force, and induction reduce to the forms 

V-(pu)=0, u ■ Vu + + n X (f2 X r) + -VP = 0, and B = /tpu, (4) 

P 

where the parameter k is constant along streamlines. Note that the magnetic fleld does not 
contribute a force along the fleld lines and that the Coriolis force also does not contribute. 
The correction to the effective gravity due to the rotating reference frame is included. 
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The velocity vector u follows the magnetic field fines. In the following sections, we 
construct coordinate systems such that one of the coordinates (denoted here as p) follows 
the field lines. As a result, the flow velocity has only one component, which points in the 
direction of the magnetic fleld p = hpe^. We also assume that the thermodynamics of the 
flow can be modeled using a polytropic equation of state of the general form 

P = Xp^+V" , (5) 

where n is the polytropic index. As shown below, however, the index n must be relatively 
large in order for the flow to pass smoothly through the sonic transition. As a result, we 
consider the limiting case of an isothermal equation of state n ^ oo for much of this paper. 



3.1. Dimensionless Formulation 

The problem can be re-formulated using dimensionless variables. Length scales can be 
measured in terms of the stellar radius We also define a reference scale pi for the density 
and use the corresponding sound speed ai as a reference velocity, given by af = [dP/dp]i, 
or, equivalent ly, 

al^K(l + ^)py\ (6) 



We can then define the following dimensionless quantities 



p lul r ^ 

— , u= — , ^=-5-^ and V'=^ 

pi ai it* ai 



Note that u is the (sonic) Mach number for isothermal flow. The stellar radius is typically 
i?* ~ 10^^ cm. The reference sound speed ai corresponds to that of the flow, which is 
launched from the inner edge of the disk; the appropriate temperature is that just above 
the disk surface, where the flow is exposed to UV heating, so that T ~ 10^ K and ai ~ 10 
km/s. The reference density scale pi depends on the mass accretion rate M; for the standard 
value M fa 10~^Mq yr"-*^, and for typical values for the other quantities (see Section 5.3), 
the reference density is given by pi/m^ = ni ~ 3 x 10^^ cm~^. 

Next we define a dimensionless parameter b that measures the depth of the gravitational 
potential well. 



R^a\ y^.bMQj cm J VlOkms 
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The depth of the potential well is thus expected to lie in the range b ~ 500 — 1000. Finally, 
we define a dimensionless parameter u that measures the effects of rotation, 

where Q is the rotation rate of the star and Ql = GM^,/Rl is the Keplerian rotation rate 
due to stellar gravity evaluated at the stellar surface. As a result, the rotation parameter 
can be written as u = b (-R*/rco)^, where Vco is the co-rotation radius. T Tauri sta rs exhibit 



a wide range of rotation periods P = 2n/n ~ 2 — 12 days (iHerbst et al.l 120071 ). For the 
median value P ~ 6 days, the co- rotation radius falls at ~ 8i?*, so we consider the range 
''^co ~ 6 — 10-R*. 

For simplicity we generally assume that the co-rotation radius r^o is coincident with the 
inner disk edge ra (so that u = b/^^, where = ra/R* ~ 6 — 10). In practice, accretion 
from the disk surface is launched from an annulus of finite extent Ar, although the width is 
generally narrow (Ar < r; see Section 5.3). If Ar ^ r, this finite width is not an issue. If 
the annulus extends beyond the co-rotation radius, the fiow falls in the "propeller" regime 
where the star is spun down (Ustyugova et al. 2006). On the other hand, if the annulus 
extends inside the co-rotation radius, magnetic links to the regions of the disk spinning faster 
than the star, as well as accretion of disk material, act to spin-up the star, which will then 
increase its rotation rate and move the co-rotation radius further inward. However, induced 
changes in the location of the corotation radius, and/or the stellar spin, generally occur on 
long time scales. 

In terms of the dimensionless fields defined above, the continuity equation takes the 

form 



du da au d 
dp dp hqh^ dp 
and the force equation becomes 



dp a dp dp dt, dp t, dp 

These equations can be integrated immediately to obtain the solutions 

auhqh^ = A , (12) 

and 

-u^ + na^^"" + ^ = e + uj [ ^sme(x-p) =e + ujl , (13) 
2 J \Vp\ 

where the second equality defines the integral / (this integral can be written in simple form; 
see Result 2 of Appendix [A]). The quantity hghip in equation ( fT2l) is proportional to the 



- 12 - 



inverse of the magnetic field strength (consistent with the third part of equation [3]). The 
parameters A and e are constant along streamlines, but are not, in general, the same for all 
streamlines (they are functions of q). 



3.2. Transitions through the Critical Points 

Critical points in the flow arise when the fluid speed is equal to the transport speed. In 
general, magnetic media support three types of MHD waves and hence allow for three types 
of critical points. In the case where the flow is confined to follow magnetic field lines, the 
problem has only one possible critical point, which occurs where the flow speed equals the 
sound speed. In order for the flow to pass smoothly through this sonic point, only particular 
values of the constant A are allowed. The equations of motion (fTOj) and (fTT]) provide two 
equations for the two unknowns du/dp and da /dp] by solving for these quantities and 
requiring that the functions are continuous at the sonic point (e.g., Shu 1992), we find the 
required matching conditions 



2 = ^1/'^ and ^:^[h^h^)+uisiTid\Vp\-^{x-p) = (14) 



u 



hqhfj, dp ^"^ dp 



Note that these expressions remain valid in the isothermal limit where n — )■ oo and hence 
^i/n _s, -|^_ must thus evaluate the geometrical factor Q defined by 

^ hqh(f, dp^ hp dp ' ^ ^ 

where the second equality follows from Result 1 in Appendix |X1 Note that for radial flow 
in spherical coordinates, where p = C,, this factor would have the familiar form Q = 2/^. 

In general, the divergence operator, and the geometrical factor Q have units of p^^, but 
the coordinate p is not necessarily a length scale. Nonetheless, we can define a dimensionless 
parameter T that measures the degree of divergence through the construction 

T = i:{u) = ^H (^^y g ^(0 = ^^' (16) 

where the function H is defined by equation f lA16|) in Appendix |Al The intermediate ex- 
pression is general, whereas the final equality applies to the specific case where the flow is 
confined to a field line so that the angle 9{^) is a known function of S,- With the defini- 
tion f ll6p . the index T is of order unity and is a slowly varying function of the coordinates. 
Similarly, we define the function A that specifies the rotational term in the force equation, 

K = ^sme[^^ (x-Vp). (17) 
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In terms of these functions, the matching condition at the sonic point takes the seemingly 
simple form 

ay-l + coA = ^. (18) 

In general the quantities T and A are functions of 6). For flow along a magnetic field line, 
which are lines of constant g, the angle ^ is a specified function of the radial coordinate, so 
that we obtain functions of a single variable T(,^) and A(^). In the isothermal limit n — t- oo, 
the density dependence drops out, so that all of the terms in equation ( ITSl) become known 
functions of the variable ^. For a general polytropic index, the two equations of motion 
(evaluated at the sonic point), in conjunction with the two matching conditions, allow one 
to solve for the four unknowns ^s-, ots-, Us, and A. After eliminating the other three variables, 
the equation that specifies the sonic point has the form 



This equation applies to steady polytropic flow that follows any magnetic field configuration, 
where the geometry is encapsulated in the index T(^) of the divergence operator, the scale 
factor hp, the rotational term A(,^), and its integral /. The physical system is specified by the 
depth b of the potential well, the rotational speed (through u), and the polytropic index n. 
Finally, the streamline of the flow is determined by the intersection point with the equatorial 
plane, where this point (often taken to be the inner disk edge ^d) specifles the value of the 
coordinate q. Equation (fT9|) remains valid in the isothermal limit n — )■ oo, where the right 
hand side of the equation becomes unity. 



A = ftj'!^l , (20) 



With the location of the sonic point specifled by equation (1191) . the dimensionless 
mass accretion rate A for flnite n is given by 

n+1/2 

where all of the quantities are evaluated at C,s- In the isothermal limit, accretion parameter 
A is specified implicitly through the transcendental equation 

\n\-^\^ = \nhp + ujl + b(^^~?j-^. (21) 



3.3. General Constraint on Steady Polytropic Flow 



In this section, we derive a general constraint that must be met for steady transonic 
accretion to reach free-fall speeds in the inner limit. Here, the flow starts at the inner disk 
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edge (at subsonic speed) and ends on the stellar surface (at supersonic speed). This result 
applies to polytropic flow that follows magnetic field lines. Let the integer £ denote the 
highest order multipole of the field near the stellar surface and let n denote the polytropic 
index of the equation of state. We find that self-consistent transonic solutions must satisfy 
the requirement 

n>£+^ (22) 

in order to achieve (nearly) free-fall speeds in the inner limit. This constraint implies that, 
to leading order, the fiow can be described using an isothermal equation of state. 

The constraint of equation fl22]) can be derived as follows: In the limit ^ <^ 1, only the 
highest order multipole field contributes to the magnetic field. In this limit, the continuity 
equation reduces to the form 

au = AAr^'+'^ , (23) 

where A is a constant (for example, this constant has values A = 2, 2r, and 2 for the 
magnetic field configurations considered in Sections H] and [5l and Appendix [Bl respectively). 
The force equation reduces to the form 

+ na'/^ = \ . (24) 
2 4 

The combination of these two reduced equations of motion results in the expression 

-u" + n (AA)^/"r^''-'^/"w"'/" = \ ■ (25) 
2 

In order for the solutions to approach a free-fall form in the inner limit, the velocity field 
must approach the form v? ~ 26/^, which requires the first term to dominate the second in 
equation (|25l) . Consistency thus requires 

^2 > a^l^ which imphes ^ > -^(i^j:^ ■ (26) 

The requirement implies that £ + 3/2 < ra, as claimed in the constraint of equation 0221) . 

We note that this argument applies for the limit where ^ — )■ 0. In practice, magnetically 
controlled accretion takes place over a limited range in radius, from the inner disk edge to 
the stellar surface, i.e., spanning only about a factor of 10 in radial scale. In order for the 
equation of state to be stiff enough to affect the accretion flow over this more limited range, 
the polytropic index n must be larger than the value £ -|- 3/2, as given by equation (12^ . 

In addition, for the case of dipole flelds = 1, Koldoba et al. (2002) flnd that the behav- 
ior of accretions flows from disks onto magnetized stars has qualitatively different behavior 
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when the polytropic index n is smaller or larger than n = 5/2 = £ + 3/2 (the same threshold 
indicated by equation [22]). For n > 5/2, the Mach number increases monotonically as 
the flow approaches the stellar surface (i.e., for decreasing radius); for n < 5/2, the Mach 
number first increases and then decreases as the fiow moves inward (see Koldoba et al. 2002 
for further discussion). 

The observational implications of this constraint are important: First we note that cur- 
rent observations indicate that transonic flow with free-fall speeds does take place. Further, 
observations of magnetic fleld signatures for T Tauri stars suggest that the octupole compo- 
nent provides the largest contribution near the stellar surface. As a result, we should take £ 
= 3, so that the constraint becomes n > 9/2. From the inner disk edge to the stellar surface, 
the density can vary significantly, e.g., by a factor of ~ 100 for the models considered in 
this paper. Due to this constraint, however, the temperature (sound speed) could vary by 
at most a factor of ~ 2.8 (1.7). As a result, the fiow is close to isothermal, and we consider 
the isothermal limit {n — )■ oo) for most of this w ork. This approximation is consistent with 
previous modeling work (IHartmann et al.lll994j ). which finds that reasonable hydrogen line 
profiles can be obtained if the temperature profiles are roughly isothermal (except for an 
initial rise in temperature where the fiow leaves the disk). On a related note, calculations of 
the thermal structure of accretion funnels (Martin 1996) indicate that the Ca II and Mg II 
ions have a strong cooling effect on the fiow, i.e., they act as a thermostat and give rise to 
nearly isothermal conditions in the vicinity of the stellar surface. 



3.4. Mass Accretion Rate 

The discussion thus far has focused on one streamline at a time. This section shows how 
the streamlines add up to determine the mass accretion rate from the disk onto the star. 

The total mass accretion rate leaving the inner portion of the disk takes the form 

/•r2 /•!;2 r£,2 

Md = 2 / 2'Krdr pv = AnRlpiai / zuduj aiUi = AnRlpiai / zudzuX (hp^)^, (27) 

Jra J(d J£,d 

where is the inner disk edge, the accreting annulus extends from to r2, and the variable 
w = r / R^. The leading factor of 2 arises because material accretes from both the top and 
bottom surfaces of the disk. In the final equality we have used the boundary conditions 
that a = 1 and u = X/hp at the start of the trajectory, and the subscript indicates that the 
quantity is to be evaluated in the disk plane. The fiow is assumed to leave the disk surface 
in the vertical (z) direction, consistent with the magnetic field line geometries of interest 
(primarily dipoles and octupoles). For the sake of completeness, however, we note that the 
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flow is somewhat more complicated because the magnetic field threading the disk bends 
away from the vertical at the disk surface, and because the surface itself is not perfectly 
fiat. As a result, material in the accretion fiow must first climb over a potential maximum 
before leaving the disk and falling toward the star (Scharlemann 1978, Ghosh & Lamb 
1979). In practice, of course, we see ample observational evidence for accretion columns and 
magnetically truncated disks (see Section [1] for references), so that material can readily leave 
the plane of the disk. 

The streamlines, which follow the magnetic field lines, end on the stellar surface, where 
the mass accretion rate M* onto the star is given by 



where the subscripts denote that the quantities are to be evaluated at the stellar surface. The 
end points of integration fii and /i2 (where n = cos 6) are determined by the angles where 
the streamlines (field lines) starting at and ^2 intersect the stellar surface. The leading 
factor of 2 arises because accretion onto the star takes place into rings on both hemispheres. 
Notice that the dot product (f ■ p) must be included because, in general, the streamlines do 
not intersect the stellar surface in the radial direction. 

In order for mass to be conserved, as required, the mass accretion fiow leaving the disk 
surface, given by equation fl271) . must be the same as the mass accretion rate onto the stellar 
surface, given by equation ( l28|) . Comparison of the two expressions shows that they are 
equal if the integrals are equal. This equality holds and is shown in Appendix O (which thus 
provides a consistency check on this approach). 



This section considers the field geometry to be that of a simple dipole. Although mag- 
netic accretion with a dipole field has been addressed previously (e.g., Hartmann et al. 
1994, Koldoba et al. 2002), we revisit the problem using our approach, which sets up the 
generalizations in the following section. 




(28) 



4. Magnetic Accretion with a Dipole Field 



4.1. The Coordinate System 



For the case of a pure dipole field, we use coordinates of the form 



p = ^ cos 9 and Q = ^ ^ sin^ 9 



(29) 
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where ^ = r/i?* (see also Radoski 1967, Canalle et al. 2005, and Appendix The 
streamline (and hence magnetic field line) that intersects the disk at its inner edge is thus 
given by g = i?*/rrf, so that this streamline corresponds to the trajectory ^-R* = r^sin^^. 

The covariant basis vectors ej arise from the gradients of the scalar fields that define 
the coordinates. If we express these vectors in terms of the original spherical coordinates 
(^, 9, 0), the basis takes the form 

ep = 2C^cosef + C^sm9e, (30) 



and 



-9 

1 



= -r^ sin^ er + 2^^ sin 9cos99, (31) 

, ■ nh (32) 

^ 4 sm 

where the quantities (f, 9, (p) are the usual unit vectors for spherical coordinates. The scale 
factors thus become 

hp = [4 cos^ 9 + sin^ 9] '^''^ , (33) 

K = [4 cos^ 9 + sin2 9] '^^^ , (34) 
smfe' 

and 

hff, = sin 9 . (35) 
The geometrical factor Q in the divergence operator is thus 

^9 ^8 — 5 sin^ 9 , . 

g = 3ecos9 . (36) 

(4 — 3 sm 9Y 

For fiow along a field line (constant g), the effective index T of the divergence operator (see 
equation [TB]) thus takes the form 

T = T(0 = |f|^V (37) 



2 V4-3«? 

In the inner limit where ^ — )■ 0, the index T — 3, as expected for dipole geometry. For 
streamlines, and hence magnetic field lines, that connect the disk plane to the stellar surface, 
the product — )• 1 where the streamlines intersect the disk. As a result, at the outer starting 
point of a streamline, the index approaches the not-so-obvious value T — )■ 9/2. Finally, for 
dipole field configurations, the rotational term A (see equation [17]) takes the simple form 



A = m = 



(38) 
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4.2. Isothermal Accretion Flow 

We now consider isothermal accretion flow. For a given streamline labeled by the co- 
ordinate g, we can eliminate the angular dependence from the equation of motion. Even 
though the fluid flelds are functions of the coordinate p only, the resulting expressions are 
simpler in terms of the radial variable ^{p)- This change of variables is transparent as long 
as ^ is monotonic as a function of p. The equations of motion thus take the form 

du da T 3 8 — 5q^ . . 

+ u—- = —au — = —au — -, (39) 

and 

du 1 da h . 6 3 , 
The matching condition for the sonic point takes the form 

3^8-5ge ^26-3a;ge'- (41) 
^ 4 - 3gC ^ ^ ^ ' 

For the case where the inner disk edge corresponds to the co-rotation point, the solution to 
equation f HTj) lies just inside the point = (2/3)^^'^ fd/ R*. We denote the sonic point as ^s- 

The equations of motion can be integrated to take the forms 

au = \{qi)-\A-?,qif^ , (42) 

and 

-M2 + lna = £ + - + -wg^^ (43) 

At the inner disk edge, ^ = 1/q and a = 1; the integration constant A in equation (H2|) is 
defined so that u = A at this boundary. The second constant e from equation ( H3ll is then 
given by 

e = h'-hq, (44) 

where we have assumed that the inner disk edge, the launching point of the flow, and the 
CO- rotation radius coincide (so that u = bq^). Using this result to specify e and taking the 
logarithm of the continuity equation fj42|) . we flnd the following implicit specification of the 
remaining constant A: 

InA - = 31ne. - ^ ln(4 - 3^.) " ^ + ^ + - ^ • (45) 

This equation has two roots for A, one with A < 1 and another with A > 1. The accretion 
solutions of interest here start with subsonic speeds at ^ = C,d = ^/q, where u = A, so we 
must take the smaller root (A < 1). 
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Fig. 2. — Dimensionless fluid fields for accretion solution using a pure dipole geometry and 
an isothermal equation of state. The inner disk edge lies at ^ = = and the stellar 
surface corresponds to ^ = ^* = 1.0. For dimensionless depth of the gravitational potential 
well b — 500, the velocity field u is shown by the solid curve and the density a is shown by 
the dashed curve. The corresponding solutions for b = 1000 are shown by the dotted curves. 
For the reference values rii = 3 x 10^^ g cm""^ and ai = 10 km s~^, these solutions have 
physical speeds -u* ^ 290 km s~^ (410 km s^^) and number densities n* ~ 2 x 10^^ g cm~^ 
(5 X 10^^ g cm"^) at the stellar surface for b = 500 (1000). 
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For flow that passes smoothly through the sonic point, we can find the dimensionless 
fluid fields a and u. These profiles are plotted in Figure |2] for two choices of the dimensionless 
depth b of the gravitational potential well of the star. Notice that the density field initially 
decreases inward from the starting point (at the disk edge), but eventually increases. This 
behavior can be understood by finding the limiting form of the solutions for the regimes 
|1 — ^1 <^ 1 and ^ ^ 1, as shown in the following subsection. 



4.3. Limiting Forms for the Flow Solutions 



In the regime |1 — g^l <^ 1, we define r] such that qC, = l — r], where ?7 -C 1 (and q = l/C,d) 
The equations of motion take the form 



and 



au = X { 1 + ^1] 



^■u^ + In a = (A^ - bqr]) 



(46) 



(47) 



where we have kept only the leading order terms (in rj). Working to leading order (in rj), we 
find 



u ^ X exp 



-ri{bq + 9) 



and 



a 



1 9 

1 + -r/ ) exp 



(45 



These expressions can be expanded further, keeping only the leading order terms in t], to 
obtain the forms 



M ~ A 



1 + ^^{bq + 9) 



and a ~ 1 — -bqr] , 



(49) 



The initial decrease in the density field is thus clear (see Figure |2]). 

In the opposite limit where ^ ^ 1, the equation of motion reduce to the forms 

au = 2Xq'^C^ , (50) 

and 



-u + In a = - . 

To leading order, the dimensionless fluid flelds become 



u 



2b 



1/2 



and a ~ 



1/2 



Ag-3 

^5/2 



(51) 



(52) 
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The fluid fields thus increase as — 0. 

One can also show that the velocity field u{^) is monotonic for this problem. If we use 
the equations of motion to eliminate the density a in favor of the velocity u, we obtain 

-u"^ -\nu = e + ^ + -bq^^^ -\nX + 3\n{q^) --\n[4:-3q^] . (53) 
2 42 2 

Taking the derivative of both sides we find 

/ 1\ du b 3 , 4 ^2 3 3q , ^, 

One can show that the function F{C,), the right hand side of the above expression, is mono- 
tonically increasing on the interval < < 1, has a zero at the sonic point, is positive for 
^ = 1 and is negative in the limit — ?■ 0. These properties, in conjunction with equation 
f l5^ imply that u is monotonic and decreasing. Note that the statement that u is decreasing 
means that the velocity increases as fluid elements approach the star. 



5. Magnetic Accretion with Dipole plus Octupole Field 

In this section we consider the stellar magnetic field to have both dipole and octupole 
components. The magnetic field thus takes the form 



3 

(5 cos^ - 3) cos f + - (5 cos^ e-l)sme9 



+ -^r^ (2 cos ^ f + sin e ej 



(55) 

where ^ = r/R^ is the dimensionless radius. The leading factors of 1/2 (for both the dipole 
and octupole terms) are included to be consistent with the convention of Gregory et al. 
(2010). If we scale out the dipole field strength, the relative size of the octupole contribution 
is given by the dimensionless parameter 



Boct 



B. 



dip 



(56) 



Observations of the signatures of magnetic accretion onto T Tauri stars indicate that the 
parameter F lies in the range < F < 10. For example, the young stars V2129 Oph (Donati 
et al. 2007) and BP Tau (Donati et al. 2008) are observed to have field parameter F = 1 — 4. 
The star AA Tau (Donati et al. 2010b) has a nearly dipole field with F ^ 0.25, whereas TW 
Hya has a much larger octupole component with F 4 observed at one epoch and F 6 
found at another (Donati et al. 2011b). 
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5.1. The Coordinate System 



With the configuration of equation fl55l) . the magnetic field is current-free and curl- free, 
and can be written as the gradient of a scalar field (see Appendix . The first scalar field 
p of the coordinate system takes the form 

1 



P 



-C^T (5 cos^ - 3) cos - cos 6 . 



(57) 



The gradient Vp points in the direction of the magnetic field. Next we construct the per- 
pendicular vector field Vg, where the scalar field q provides the second coordinate and is 
given by 

-'^ (5^ 



1 



The scalar fields {p, q) represent orthogonal coordinates in the poloidal plane and are used 
instead of the spherical coordinates (^, 6). In this version of the problem, the magnetic field 
is axisymmetric about the z axis, so that the usual azimuthal coordinate is the third scalar 
field. Note that both p and q are dimensionless. 

Next we find the covariant basis vectors e^, which can be written in terms of the original 
coordinates (^, 9), so that the basis takes the form 



[r^T (5 cos^ ^ - 3) + 2r^] COS 9 1 



-r'r (5cos2 0-l)+e 



-3 



sin^^, 



(59) 



and 



jr^r (5cos^0-i)+r 



sin^ 9f+ [C^r (5 cos^ ^ - 3) + 2^^] sin 9cos9 9, (60) 

where the third basis vector is given by equation 
It is useful to define ancillary functions 

/ = r (5cos2^-3) +2^^ and ^ = (5 cos^ - l) + (61) 
With these definitions, we can write the magnitudes of the basis vectors in the forms 

k/ = r''[fcos'9 + g'sm'9] , (62) 

and 

|e^|2 = sin^ 9 [g^ sin^ 9 + f cos^ 9] . (63) 
The corresponding scale factors thus become 



= ^5 [f cos^ ^ + sin^ 



-1/2 



(64) 
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and 

K = ^ [9' sin^ e + f cos' e] , (65) 

where the third scale factor is given by equation ( l35ll . 

This specification of the divergence operator is imphcit. One could invert equations 
f l57|) and f l58|) . and then write the spherical coordinates (^,6^), the scale factors {hp, hg, h^), 
and the ancillary functions {f,g) as functions of the new coordinates {p,q)- However, the 
definitions of {p,q) are nontrivial functions of {C,,0), so that inversion is complicated and 
unwieldy. For clarity, we leave this construction in implicit form. 



5.2. Results for Flow Along a Field Line 

The lines of constant q define the trajectories of fiuid elements following along magnetic 
field lines. Here we are interested in the subset of field lines that intersect the equatorial 
plane and hence intersect the disk. Note that in a fixed poloidal plane, and close to the 
stellar surface, the axisymmetric dipole plus octupole fields have three regions of closed field 
line loops between the north and south pole of the star; the higher latitude regions that do 
not intersect the disk are not considered here. By inverting equation (!58|) . we can find the 
angle 9 function of radius ^ along a field line: 

sin^ = ^{{e + r)- [ie + rr - 5rqe] '^'} . (66) 

For the streamlines of interest, we must choose the negative sign for the discriminant (as 
shown below in equation [7T], for sufficiently large values of F, streamlines that intersect the 
disk have q < 0, and we must use the opposite sign). The ancillary functions can then be 
written as a function of the radial coordinate only: 

f(^) = 2[ie+rr-5TqeY^' (6?) 

and 

9i0 = ^ {r + [ie + - 5rqe] '^'} - . (68) 

With these definitions, we can find the rotational function A (see equation [T7]) for fiow 
along a field line, 

m = § {(e' + n - [ie + rr - ^vqey^'} {i + f§} • (69) 

The index T (from equation [I6]) takes the form 

-1 

(70) 



T(0 = 5 - 



o 
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Fig. 3. — Magnetic field lines centered on the star for configurations with both a dipole and 
an octupole contribution. For this case, the dimensionless parameter F = 10. The star is 
shown in the lower left of the figure, the disk is shown as a wedge just above the x-axis. A 
collection of hues of constant p, the surface orthogonal to the field hues, is also shown. 
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5T 



9 [^h-^] - fh 



{2e +2T-f) + - n (2^ - i/,) 



where /(^) and g{^) are given by equations fj67t [681) . and where = df /d^. 

The index T is a slowly varying function of the radius ^, or, equivalently, the position 
p{S,) along the field line. The index function T is plotted in Figure HI for varying values 
of the parameter F that sets the relative strength of the octupole component. The solid 
curves show the index T(^) for parameter values F from 10~^ (bottom) to 10'^/^ (top). The 
dashed curve shows the limiting case of a pure dipole field. In the limit ^ — )• 0, the index 
T — )■ 5, as expected for an octupole field. For a pure dipole field (dashed curve), the limiting 
value T — 3. Near the inner disk edge, the index is somewhat larger than the value T 
= 9/2 expected for a dipole. For the intermediate regime in ^, the index first decreases 
as decreases (toward the dipole value of 3), but the index increases closer to the star as 
the octupole contribution dominates. For sufficiently large values of F ^ 1, the octupole 
component dominates the field even near the inner disk edge; for this regime, the value of 
the index approaches T = 15/2 (from equation [70] in the limit F — > 00). Notice also that 
for the value F ~ 30 (e.g., the top curve in Figure H]), the compromise between the dipole 
and octupole terms leads to the index being nearly constant with value T ^ 5. 

The field lines cross the equatorial plane where sin 9 = 1 and cos ^ = (where the fiuid 
trajectories start). The value of q that leads to plane-crossing at radius is thus given by 



9 = - ( 1 - 772 ) ■ (71) 



The dimensionless truncation radius ~ 10, so we expect q ~ 1/10 for field lines that 
thread the disk. If the octupole contribution is too large, F > ~ 100, the field lines that 
intersect the disk would have g < and the sign of the discriminant in equation 0661) must 
change. T Tauri systems are expected to have F < 10, so this complication should not arise. 

Using the result g ^ 1, we can find an approximate expression for the co-latitude on 
the star where fields lines starting at the disk truncation point reach the stellar surface, i.e.. 

The second equality uses equation (I7T]) to write q in terms of the disk truncation radius 
^d- This result shows that the field lines that connect the disk truncation to the star must 
meet the stellar surface near the pole (^* ^ 1). Increasing the octupole contribution (via 
increasing F) decreases the angle 9,, even further. 

We can also find the ratio of areas. Let Ad be the area of an annulus on the disk. In 
the limit of a thin annulus with width A,^, the area Ad = 27i^{A^). This area gets funneled 



- 26 - 




Fig. 4. — Index T of the divergence operator (see equations [I6] and [70]) for magnetic 
field configurations with both dipole and octupole components. The relative strength of the 
octupole component T has values given by logj^gF = 3/2, 1, 1/2, 0, -1/2, and -1 (from top 
to bottom). The dashed curve shows the limiting case of a pure dipole field where F 0. 
The heavy solid line at ^ = 1 marks the surface of the star, where accretion flow stops; the 
curves continue inward to smaller values of ^ to illustrate the asymptotic behavior of T(^). 
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onto the stellar surface over a much smaller area A* given by 

A* = 2tt [cos^^*2 — cos^^*i] . (73) 

This expression does not include the dot product f ■ p, which takes into account the non- 
radial direction of the incoming material (equation [28] and Appendix [U]). Since the angle 
9^ <^ 1, the dot product is close to unity, r ■ p ^ 1 — (SF + 1)^/ (8.^rf[r + 1]^) + . . ., and can 
be neglected to leading order. Using the above results, we can evaluate the ratio of areas to 
find 

This leading order expression shows that the ratio of areas decreases inversely with the 
increase in field strength, which increases as from the dipole contribution and includes 
an extra factor of (1 + F) at the stellar surface due to the octupole component. 

In the limit of small ^, we can find asymptotic forms for the velocity and density fields. 
As long as the equation of state is not too stiff, the dimensionless fluid fields approach the 
form 

u ~ (^j^ and « ~ 0) A (r + 1) r'/' • (75) 

These forms are valid when the polytropic index n > 9/2 (see Section [33]) . For flow with 
n < 9/2, free-fall velocities are not realized. 



5.3. Accretion Solutions 

Figure shows the density and velocity fields for a typical system with magnetic oc- 
tupole parameter F = 10 and dimensionless depth of the gravitational potential b = 500. 
The solid curve shows the dimensionless velocity field u{^) and the dashed curve shows the 
corresponding density profile a(^). For comparison, the dotted curves show the profiles for 
the same system with a pure dipole field. Since the dipole contribution dominates the mag- 
netic field at the inner disk edge, and since the sonic point falls relatively near the edge, the 
location of the sonic point and the mass accretion constant A are only weakly dependent 
on the parameter F. Specifically, the sonic point C,s ~ 8.80 for F = 10, compared with C,s ~ 
8.87 for F — )■ 0; similarly, the parameter A ~ 0.0877 for F = 10, compared with A ~ 0.113 for 
the limit F — )■ 0. Notice, however, this statement no longer holds if the octupole component 
dominates the magnetic field at the inner disk edge. 

The most important difference between the two field configurations (with solutions 
shown in Figure [5]) is that the density increases near the stellar surface more rapidly in 
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Fig. 5. — Dimensionless fluid flelds for the velocity u (solid curve) and density a (dashed 
curve) for accretion solution with both dipole and octupole components. In these units, 
the inner disk edge corresponds to ,^ = = 10 and the stellar surface lies at = = 1. 
The relative strength of the octupole component is set by the parameter F = 10, and the 
dimensionless depth of the gravitational potential well h = 500. The dotted curves show the 
corresponding solution for a dipole field (in the limit F — )■ 0). For the reference values rii = 
3 X 10^^ g cm~^ and oi = 10 km s~^, these solutions have speeds 290 km s~^ at the 
stellar surface; the corresponding number densities are n* ~ 2 x 10^^ g cm""^ for the pure 
dipole and 2 x 10^^ g cm~^ for the dipole/octupole configuration. 
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the presence of an octupole component; here the density is larger than that of the dipole 
hmit by a factor of ~ 9.5. Equation fl75|) suggests that the density should be larger than that 
of the dipole case by a factor of (1 + F) = 11; this factor is somewhat larger than the value 
(9.5) depicted in Figure E] because the solutions are not fully in the ^ — )■ limit. Nonetheless, 
the inclusion of octupole field components allows the flow densities at the stellar surface to 
be greater by an order of magnitude compared to those from the dipole limit. 

The physical value of the flow density at the stellar surface, before the accretion shock, 
can be estimated as follows. The mass accretion rate can be written in the approximate 
form 

M = 27rpiai(r2 + rrf)(r2-rd)(A), (76) 

where (A) is normalized so that m = A at the inner disk edge and the mean value is taken 
to account for possible variations over the range of launching radii on the disk surface. 
It is useful to define a dimensionless width of the annulus where the flow originates, i.e., 
w = {r2 — rd)/rcL- The number density ni at the start of the flow can then be written in the 
form 



3 X 10^^ gem ^ 



M 



[l + w/2)w \lO-8M0yr 



-1 



ai 



10 kms 



-1 



-1 



(A) 
0.1 



10^^ cm 



(77) 



For dipole field configurations (see Figure |2]), the number density at the stellar surface 
is larger than the initial density rii by a factor of 1 — 10, so that n* ~ 10^^ g cm~^ (where 
we have taken (1 + w/2)w ~ 1). For field configurations with octupole contributions, the 
number density is larger by another factor of ~ (l + F) (see Figure E]), so we expect ^ 10^^ 
g cm~'^ for systems with substantial octupole components. The densities could be even larger 
if the initial annulus on the disk is narrow (w -C 1). 

For comparis on, coronal densi t ies for T Tauri star s are estimated to be of order ricor ~ 



10^ — 10^^ g cm ^ (INess et al.ll2004l : iKastner et al.ll2004l ). and these values are considered too 



low t o produce the observed soft X-ray exess emission from these sources (IBrickhouse et al. 



3x10^ 



20101 ) . The X-ray emission can arise from sh ock heated plasma wi th temperatures T 
K and densities n ~ 10^^ — 10^^ g cm~^ (lArgirofE et al.ll201ll ). These required densities 
are thus comparable to those expected from accretion flow (see above); more significantly, 
octupole contributions F 7^ may be necessary to explain the upper end of this range. 

The expected accretion hot spot temperature T5 can be written in terms of the other 
system parameters such that 



GM^M 



1 

Td 



2 ^J(i + r) gmM 

Tc w{2 + w) asR^rj 
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where is the Stefan-Boltzmann constant, and the other quantities have been defined 
previously. Note that T5 is the hot spot temperature at the end of the accretion flow where 
material arrives at the star and generates an optically thick shock. For the same physical 
values used to evaluate equation (1771) . we obtain T5 ^ 7000 K (1 + F)^/^. Even for a 
strong octupole component, F ~ 10, the hot spot temperature is only Tg ~ 12, 700 K. As 
a result, the hot spot temperatures for accretion flow along non-dipolar magnetic field lines 
remain consistent with the constraints established by Muzerolle et al. (2001) from line profile 
modeling, line ratios, and continuum emission constraints (see their Figure 16 and associated 
discussion). We note that the higher temperatures derived form X-ray observations of He- 
like line triplets (e.g., Argiroffi et al. 2011), and discussed in the previous paragraph, refer 
to the hotter post-shock regions. Although a full treatment of X-ray signatures is beyond 
the scope of this present work, other studies have suggested that emitted X-rays from the 
accretion shock could be absorbed within the accretion funnel (Sacco et al. 2010). 

Figure [6] shows the mass accretion parameter A as a function of the depth h of the stellar 
gravitational potential well. Curves are shown for a range of values for the parameter F that 
sets the relative strength of the octupole component of the magnetic field. As shown in the 
figure, increasing the strength of the octupole contribution results in a decrease in the mass 
accretion A. For all values of F, the parameter A decreases in both the limit of large h and 
the limit of small h. For sufficiently large 6, Figure [6] shows that the accretion parameter 
A is an exponentially decreasing function of the gravitational potential h. This behavior 
can be understood as follows: In this regime, the matching condition for the sonic point 
from equation f[T^ reduces to the form uKS,'^ = b; since u oc b, the location of the sonic 
point becomes independent of the potential b. The mass accretion parameter is specified 
through equation (12T1) . which takes the form \nX = A — Cb, where A and C are constant 
and C = 1 — 1/^5 > 0. The mass accretion parameter is then given by A ~ exp[—Cb], as 
depicted in Figured 



6. Magnetic Truncation Radius 

In these accreting systems, the circumstellar disk is (indirectly) observed to have an inner 
boundary at radius > i?*. Further, the inner disk can be truncated at radius tt where 
the stress due to the stellar magnetic field is large enough to remove angular momentum 
from the Keplerian fiow (for further discussion, see Ghosh & Lamb 1979, Konigl 1991, Shu 
et al. 1994). Here we assume that = vt, but more complicated possibilities remain. The 
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Fig. 6. — Dimensionless mass accretion parameter A as a function of the dimensionless depth 
b of the gravitational potential well of the star. The magnetic field has both dipole and 
octupole contributions, where the parameter F specifies the relative strength of the octupole 
component. The curves show X{b) for varying values F = 30 (dot-dashed), 10 (solid), 3 
(dotted) , and in the pure dipole hmit F — )■ (dashed) . 
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requirement of disk truncation implies a constraint on the magnetic field of the form 

1/2 



(79) 



where a is a dimensionless parameter of order unity; for example, recent numerical simula- 
tions (Long et al. 2005) suggest that a = 1/2 (see also the review of Bouvier et al. 2007, and 
references therein) . This parameter incorporates the difference between exact equality in the 
two sides, the difference between the flow speed and the free-fall speed, and the departure 
of the geometry from spherical symmetry. 

For a dipole field, we can write B{r) = (1/2) B^ip {R*/r)^ for the magnetic field strength 
at the equatorial plane. Equation fl7^ then reduces to the form 

rr = ^-2/7 (^^^^12^1 . (80) 



For typical parameters of T Tauri stars, the truncation radius = / R* = ^ ^ 10. 

For a magnetic field with both dipole and octupole components, the truncation radius 
is given by the solution to the equation 



1 - -rr' 

4 ^ 



f R* \ _ ^7/2 

" 4M [gmJ 



where the quantity ^to is the dimensionless truncation radius for a pure dipole field. For 
the standard (aligned) field orientation with F > 0, equation (|8T|) shows that the magnetic 
truncation radius is smaller for field configurations that include octupole components, even 
though the surface field strength is larger (see also Gregory et al. 2008). The truncation 
radius is thus given by 



TO 





4/7 




4 ^ 







(82) 



where ^^o is given by equation ( IHOjl . and where the second (approximate) equality is correct 
only to leading order. 

It is straightforward to see why the disk truncation radius becomes smaller in the pres- 
ence of an aligned octupole component: For a pure dipole field, in a fixed poloidal plane, 
consider a field line that connects the positive field region in the northern hemisphere to the 
negative field region in the southern hemisphere; the field vector B points downwards in the 
equatorial plane of the star (the mid-plane of the disk). For an octupole field, however, the 
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field line that crosses the equator connects a negative field region in the northern hemisphere 
to a positive field region in the southern hemisphere, so the field vector B points upwards in 
the equatorial plane (opposite to that of the dipole field). For composite field configurations, 
when the field vectors representing the dipole and octupole parts of the field are added in 
the mid-plane, they are anti-parallel, and the resultant field vector is less than that of a pure 
dipole magnetic field. Since the field strength in the mid-plane is smaller, the disk is able to 
push closer to the star. 

Notice that the octupole component can be anti-parallel to the dipole. In this case, the 
results of this paper remain valid, with the parameter F — )■ — F. Indeed, magnetic maps 
of the young star TW Hya (Donati et al. 2011b) indicate that the star has an octupole 
component significantly larger than the dipole component, and that the two contributions 
are anti-parallel (F fa —6). More specifically, the negative pole of the octupole coincides with 
the positive pole of the dipole and the visible rotation pole of the star (both components are 
roughly aligned with respect to the rotation axis). For such a configuration, the field vectors 
of the dipole and octupole components of the field are oriented in the same direction in the 
stellar mid-plane. The field strength of the two components thus add at the inner disk edge, 
so that the disk is truncated farther from the star than it would be with a pure dipole field. 

Finally, we note that — except in extreme cases — the dipole component of a multipolar 
stellar field is dominant in determining the disk truncation radius. The corrections due to 
the octupole field, as considered here, are usually relatively small, of order O (TRl/r^). The 
parameter F is often observed to lie in the range F = 1 — 10; however, a stellar magnetic 
field would require F ~ 100 in order for the octupole component to control disk truncation. 
A quadrupole field component (not considered here) could control truncation at a smaller 
relative field strength. Notice also that for extremely large accretion rates, the truncation 
radius can be much closer to the stellar radius, where the higher order multipoles are impor- 
tant. Finally, we stress that although the dipole term tends to control ry, the higher order 
components play a large, even dominant, role in guiding the accretion flow. 

7. Conclusion 

This paper has re-examined the problem of magnetically controlled accretion onto newly 
formed (or forming) stars. Our main results include the following: 

[1] We have constructed orthogonal coordinate systems {p, q) in the poloidal plane. One 
coordinate follows the magnetic held lines and hence the accretion flow. Speciflcally, the 
variable p measures the distance along the fleld lines, which correspond to curves of constant 
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coordinate q. We also construct the scale factors, and hence the differential operators, for 
these coordinate systems. This paper constructs coordinate systems for magnetic field con- 
figurations with a pure dipole component (Section Hj), for dipole and octupole contributions 
(Section E]), and for dipole and radial (split monopole) contributions (Appendix [Bl). The 
coordinates for a general multipole field can be written in terms of Legendre polynomials, 
as shown in Result 6 of Appendix |A] (see also Gregory 2011, Gregory et al. 2010). This 
approach can thus be generalized further to incorporate ever more complex magnetic field 
configurations. 

[2] Steady-state, transonic accretion cannot reach free-fall velocities if the effective equa- 
tion of state is too stiff. For dipole magnetic fields, accretion flows that appoach free-fall 
require polytropic index n > 5/2 (see also Koldoba et al. 2002). For a magnetic field geom- 
etry with higher order multipoles, with highest order given by £, this constraint is tighter 
and takes the general form n > i + 3/2 (Section 13.31) . In particular, for the octupole fields 
{£ = 3) that are inferred for many observed T Tauri star/disk systems, the index n > 9/2. 
For cases that allow transonic free-fall flow, the pressure term becomes negligible compared 
to the kinetic term in the force equation. This behavior is analogous to that found in the 
inner limit of the generalized infall-collapse flows that form the star /disk systems themselves 
(see the Appendix of Fatuzzo et al. 2004). This constraint on the polytropic index n places 
a corresponding constraint on the allowed temperatures of the flow; accreting material can- 
not increase its temperature by more than a factor of ~ 3, and hence the flow is close to 
isothermal. 

[3] This formulation of the problem allows for the location of the sonic points (when 
they exist) to be determined semi-analytically. For isothermal flow, equation (ITS!) becomes 
a function of the radial coordinate ^ only, where the index T(^) of the divergence operator 
and the rotation term A(^) depend only on the geometry of the streamlines (given here by 
the geometry of the magnetic field). The matching condition at the sonic point is given by 
equation (HTl) for isothermal flow with a dipole geometry, and by equations f|T8l [69l [70|) for 
isothermal flow with both dipole and octupole components. For this latter configuration, 
the geometry of the field is determined by the parameter F, which sets the strength of the 
octupole component of the field relative to that of the dipole. 

[4] We have used this formulation to study accretion flow with dipole magnetic field 
geometry (Section Hj) and with both dipole and octupole components (Section [5]). Compared 
with the case of a pure dipole field, the inclusion of octupole contributions funnels the flow 
onto a smaller region at higher latitudes on the stellar surface. The flow speeds and location 
of the sonic points are largely unchanged, but the flow densities are much larger (by a factor 
of ~ F) as the flow approaches the star and at the stellar surface. For systems with strong 
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octupole components, this enhancement increases the density at the stellar surface by an 
order of magnitude, roughly from n,, ^ 10^^ g cm~'^ to n^, ^ 10^^ g cm~^ (see Section 5.3). 

[5] The inclusion of higher order multipoles alters the predicted location of the magnetic 
truncation radius (Section [6]). Aligned octupole components leads to a decrease in the 
truncation radius, whereas anti-aligned octupoles increase the truncation radius. In either 
case, however, the dipole component dominates the determination of r^. 

One way to summarize this approach is by outlining the parameters of the problem: 
The accretion flow is assumed to follow the magnetic field lines, which are determined by the 
given field geometry; for the case of joint octupole/ dipole fields, for example, the geometry 
is specified through the parameter F. To find solutions to the dimensionless problem (in the 
isothermal limit), we must set the dimensionless depth h of the gravitational potential well 
and the location S,d of the inner disk edge. We thus have a three dimensional parameter 
space (r, 6, ^d)- The rotation parameter uj is specified if we assume that the corotation point 
coincides with the inner disk edge; in general, this might not hold, and u represents another 
parameter of the dimensionless problem. Conversion to physical parameters requires more 
variables to be specified: The magnetic truncation radius is a function of the stellar field 
strength B^, (in addition to F), stellar mass M*, stellar radius i?^,, and the mass accretion 
rate M. With these quantities determined, the sound speed ai is determined for a given 
value of h (see equation |8]). The width w of the accretion annulus must also be set; the 
density scale pi is then defined through equation fl77|) . 

One goal of this work was to develop coordinate systems that allow for a semi-analytic 
treatment of magnetically controlled accretion flows in complex geometries. We have ap- 
plied these results to star/disk systems with both dipole and octupole components, but this 
work should be extended in a number of directions. In terms of theoretical development, we 
have focused on the case of isothermal flow and octupole fields. The thermodynamics of the 
accretion flow can be modeled with increasing accuracy, first by using a general polytropic 
equation of state, and then by including a full treatment of heating and cooling. The field 
geometries should also be generalized, including additional multipole components and cases 
where the various magnetic poles, and the rotational pole of the star, are not aligned. This 
latter complication breaks the axial symmetry of the problem and thus requires consider- 
able development. The work presented herein is largely theoretical, so that an important 
step is to apply these techniques to specific observed sources; such modeling should also 
include comparison of line profiles. Finally, these techniques can be ap plied to a host of ad - 



ditional astrophysical problems, inc luding accretion onto w hite dwarfs (ICanalle et al.ll2005l ). 



accretion in neutrori star systems ( iGhosh fc LambI 119781) . accret i on iri black hole systems 
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trolled outflows from planets ( lAdamsll201ll ). and many others 
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A. Mathematical Results for Coordinate Systems 

This Appendix provides a collection of formal mathematical results that constrain and 
specify the class of coordinate systems used herein. 

Result 1: For the class of orthogonal coordinate systems considered here, the following 
identities must hold: 



Vg 

where 



1 dp ^ dp ^ 

— + 



\Vq? = F^\Vp\\ (Al) 



F = ^ sin ^ = /i0 , and hp = hqh^ . (A2) 

Proof: The first identity of equation f lAip follows from the requirement that Vg must be 
perpendicular to Vp. The second identity follows directly from the first. Next we show that 
the form of the function -F(^, 9) is given by the third identity (lA2p . From equation (lAip . the 
partial derivatives of q are given by 

dq Fdp dq dp 

For consistency, the partial derivatives d'^q/d^dO must be the same for either ordering, which 
implies 

d fFdp\ d f^pdp\ ^^^^ 



de\idej di V di 

If we expand and use a more compact notation, this expression becomes 

- ^ {FePe + Fpee) = F^p^ + Fp^^ + ^Fp^ . (A5) 
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The requirement of a divergence-free field, V ■ B = 0, implies that the scalar field p must 
obey the Laplace equation, V^p = 0, which requires p to satisfy the differential equation 



1 cot^ 

If we combine the previous two equations, we find 



P^^ + 7Pc + ^Pee + —prPe = . (A6) 



pe {Fe - cot OF) = ip^ {F - ^F^) . (A7) 

Since this result must hold for arbitrary field configurations, and hence for arbitrary pg and 
p^, the differential equations for F must individually vanish. These conditions imply that F 
= sin 9 as claimed (where C is a constant). The final identity hp = hqh^ follows from the 
definitions of the scale factors and the second identity. □ 

We note that this choice for F{^, 6) is not unique. One can always rescale the variables 
(e.g., q —J- Aq) or redefine the variables (e.g., q — )■ q^) and obtain a valid orthogonal coor- 
dinate system. However, the relationships between the scale factors are not invariant under 
such transformations. 

Result 2: The integral required to evaluate the rotation term can be written in the form 

/ = |esin^(x-p)^ = jmdi, (A8) 

where the integral / arises in the integrated form of the equation of motion f[T^ and where 
A(^) is defined through equation (fTTI) . 

Proof: Using the definition of A(^), we can write the integrand of / in the form 



ivpi ^^^^ ivpr 



esin^(x-p)^ = A(0 (A9) 



where we note that p = e /|e | = Vp/|Vp|. Next we write dp in the form 



86 



P,+Peg^ 



(AlO) 



dp = p^d^ + pgd9 = d$, 

The field lines correspond to lines of constant q so that 

dq = q^d^ + qgde = , (All) 

and hence 

= = (A12) 
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where the second equahty follows from the orthogonality of the coordinates p and q. Com- 
bining equations flAlOp and f lA12p allows us to write dp in the form 

|2 



dp = d^ 



Pi 



2 I 2 



d^ 



Pi 



(A13) 



Using form of the integrand from equation ( 1A9[) and the expression for dp from equation 
( IA13p . the integral I takes the form 

Np\' 



\^p\ Pi 
This confirms the result of equation (1A8I) . □ 

Result 3: The partial derivatives are given by 

pc , 89 

— = 77 and — 
H 



(AM) 



dp 



where we have defined 



Pe 

dp i'^H 
1 



h = \vp\'=pI + ±^pI 



(A15) 



(A16) 



Proof: To show the validity of this result, begin with the expressions for p(^, 9) and g(^, 9) 
and take derivatives: 

dp 89 

IT = ^=Piir +Pe^ ^ (A17) 
op op op 



and 



dq 



09 



= q^ + 
dp dp dp 



(A18) 



One can then substitute the derivatives of p for those of q (using Result 1) to obtain 

-ed-p^'% = '- ^^''^ 

Equations (]A17p and (1A19P thus provide two equations for the two unknowns d^/dp and 
d9ldp, which can be solved to obtain the stated result of equation (lAlSp . □ 

Result 4: The relationship between the derivative with respect to the coordinate p and the 
original gradient operator is given by 

d 



dp 



(A20) 



Proof: This result follows from the chain rule and the above definitions: 

d d d^ d d9 p^ d Pe d 
dp ^ didp^ d9dp^ lldi^ 'CHd9' 



(A21) 
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where we have used Result 3. This expression can be re-arranged to obtain the form 



dp 

□ 



9 



H-^/^p ■ V . (A22) 



Result 5: For any coordinate system in the poloidal plane, the components of the scalar 
fields p and q can be added: Suppose that the magnetic field has multiple components, 
represented by the vector fields Vpj, where the index j labels the component. Let qj be a 
scalar field that provides the coordinate that is orthogonal to pj. Then one can construct a 
complete orthogonal coordinate system 

Xp^J^Pj ^Q = I]?.-, (A23) 

j j 

where VXp points in the direction of the magnetic field and where VXq ■ VXp — (see also 
Backus 1988). 

Proof: By construction, we find 

v^p=i:vp,-(i:t)'-j(i:t)«. (A24) 



and 



We can write the derivatives of the qj in terms of derivatives of the pj, so that this second 
expression becomes 

VX, = Fit e)(-lj:^^y+ '-Fit e)(^iY:^-^y, (A26) 

where F{^,9) = ^sin^ is the same function for all of the components (from Result 1). 
Using this latter form for VXq, it becomes clear that VXq ■ VXp = 0. □ 

Result 6: For a magnetic field configuration that corresponds to a multipole of order i, the 
coordinates {p, q) can be written in the form 

^^_^^-(m)p^(^) and ? = 7r^sin^^^(/^)> (A27) 

where ^ is a constant, P^(//) is the Legendre polynomial of order £, and /i — cos9. 
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Proof: The form for the function p{C,, 6) follows from the requirement that the scalar field 
p must satisfy Laplace's equation and from the definitions of multipole expansions. We can 
find the form for the second scalar field 6) using either of the relations defined in Result 
1. In this context, the first of these relations takes the form 

^ = -sin^^ = -Ar(^-^)sin^^^, (A28) 

which integrates to the expression of equation flA27p . The second relation implies 

^ = esm9^ = A{e+l) sin 9 P,(/i) , (A29) 

which integrates to the same form. □ 

Note that Result 6 shows that an individual multipole component can be written in the 
form given by equation flA27p . and Result 5 shows that the field components can be added 



according to equation f lA23p . These results thus provide a blueprint to construct coordinate 



systems that trace magnetic field structures of arbitrary complexity. 



B. Coordinate System for Dipole plus Radial Field 

In this Appendix we develop a magnetic field model that includes both dipole and 
radial components, where the radial field is actually a split monopole field. For this class 
of problems, we only need to consider flow — and hence field geometry — in one quadrant 
of the poloidal plane. The distinction between a true radial field and a split monopole is 
thus unimportant for purposes of constructing the coordinate system. This type of magnetic 
field structure arises when the infall-coUapse flow that forms the disk drags in magnetic field 
lines from the original molecular cloud core. Another motivation for this structure is that 
the Solar magnetic field can be modeled using dipole, radial, and quadrupole terms (e.g., 
Banaszkiewicz et al. 1998), where the first two contributions are dominant. The coordinate 
system constructed here can thus be used to study flow along magnetic field lines for the 
Sun and other main-sequence stars with winds. 

The magnetic field under consideration takes the form 

B = B^^r^f + i^dipr' (2 cos f + sin 6 ^) , (Bl) 

where ^ is a dimensionless radius. If we scale out the dipole field strength, the relative size 
of the radial field is determined by the dimensionless parameter 

/3 = ^. (B2) 

-Ddip 
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For Solar fields, the parameter /3 must of be order unity. An even better fit to the observed 
Solar magnetic field can be obtained by using a modified split monopole term, where the 
cartesian coordinate z — )■ z + a^, and where the length scale ~ 1.5 Rq. This modification 
can be incorporated into the coordinate system, as shown below. 

For the magnetic field configuration of equation (IBip . the two scalar fields that define 
the coordinate system in the poloidal plane are given by 



(3 cos 6 Sin^ n n /-n \ 

p=-- — and q= — pcosO. (B3) 



Next we find the basis vectors 

ep = CHP^ + '^cose) r + Cosine e, (B4) 

and 

eg = -r^ sin^ ef + C^ sin 9 (2 cos 9 + (3^) 9 , (B5) 
where is the same as before. The scale factors are then given by 

hp = ^3 + 2 cos 9f + sin^ 9] '^''^ , (B6) 

and ^ 

^ = ^ [sin' ^ + (2 cos ^ + , (B7) 

where is the same as before. 

In this coordinate system, the magnetic field lines that intersect the disk must cross the 
equatorial plane at dimensionless radius = 1/g. For a given value of ^rf, the field line will 
hit the stellar surface at the co-latitude given by 



cos6'* = - 



1 r ^ / . . , . N 1/2 



-/3+(/3' + 4-4/ed)- 



fB8l 



Notice that field lines labeled by positive g > cross the mid-plane and hence are closed, 
whereas field lines with negative g < are open and extend to large distances from the star. 
On the stellar surface, the angle 9c that divides the open field lines (from the polar regions) 
and the closed field lines (from equatorial regions) is given by 2 cos 9c = 7/32 + 4 - /3. Note 
that the field lines that cross the equatorial plane make an angle A with respect to the 
vertical, where A = tan^^(/3^^); this angle must be incorporated into the analysis (e.g., the 
mass accretion rate calculation of Section 3.4). The field lines for this magnetic configuration 
are shown in Figure [71 
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Fig. 7. — Magnetic field lines centered on the star for configurations with both dipole and 
split monopole contributions. The relative strength of the monopole component is given by 
the parameter /3 = 1/4. The star is shown in the lower left of the figure. 
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Along streamlines, which correspond to lines of constant g, we can invert equation f lB3l) 
to find the angle as a function of the radial coordinate ^, i.e., 

cos9 = -^ + \[P'e + ^-^(\"'" . (B9) 

The field lines cross the equatorial plane, and intersect the disk, where cos0 = 0, which 
occurs when = 1. As a result, the field line (or streamline) that intersects the inner disk 
edge is labeled by g = 1/,^^. 



(BIO) 



Along a streamline, the magnetic field strength takes the form 

5(0 = |B| = {4 - 3ge + \p'e + \pi [P^ + 4 - 4ge] , 
and the rotation term takes the form 

m = f {2, + ^ we + 4 - 4,e] - /?^{} {3 - . (BID 

Finally, the index of the divergence operator can be written in the form 
T(0 = 3 - ^ |4 - 3ge + \p'e + \k [(i'e + 4 - 4^^] } ' x 

{-3g + + /3 [P'e + 4 - [p^e + 2 - m] } . 

(B12) 

The index T (^) is shown in Figure El 

In models of the Solar magnetic field, the split monopole contribution arises due to a 
current sheet that resides in a thin, disk-like structure (e.g., Banaszkiewicz et al. 1998). 
Outside the current sheet, we can find coordinates for the magnetic field using the methods 
of this paper. The relevant coordinates (p, q) are given by 

cos 61 M , 
P = T-o T7^, B13 

and 

^ i (e + 2eacos^ + a2)i/2 ' 
where a is a length scale that incorporates the disk-like structure and where M. sets the 
relative strength of the split monopole contribution. In the limit a — )■ 0, and the limit 
^ —7- 00, we recover the simpler split monopole model of equation ( IBip . Models of the Solar 
magnetic field imply that the length scale a ~ 1.5 (in units of the solar radius Rq) and that 
the constant Ai ~ 0.65. 
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Fig. 8. — Index T of the divergence operator (see equations [16] and |B12] ) for magnetic 
field configurations with both dipole and spht monopole components. The relative strength 
of the monopole component (3 = 2, 1, 1/2, 1/4, 1/8, 1/6, and 1/32 (from bottom to top). 
The heavy solid line at ^ = 1 marks the surface of the star, where accretion flow stops. 
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C. Conservation of Mass 

For consistency, tlie mass accretion rate leaving the disk surface and that striking the 
stellar surface must be the same. This Appendix shows that this constraint is satisfied. 
Specifically, we require that the expression for the mass accretion rate from the disk, given 
by equation fl27|) . must equal that for the mass accretion rate onto the star, given by equation 
( |28|) . The two integrals must be identical, i.e., 

z^dz^X {h;')^ = / rf/i (f.p)^ A {h;')^ , (Cl) 

where all quantities are evaluated in the disk mid-plane for the first integral and are evaluated 
at the stellar surface for the second integral. To show equality, we must find the relationship 
between the two integration variables, those on either side of equation ( ICll) . The streamlines 
are given by g = q{C,, 6), where g is a constant along a streamline. On the disk surface 6 = 
it/2, and scalar field g is a function of the radial variable only. We use w = to denote the 
radial variable for the disk surface, and A to denote the function, so that q = A{w) on the 
disk. Similarly, on the stellar surface ^ = 1, and the field g is a function of the angle only. 
Here we use the variable fi = cos 6 to denote the variable on the stellar surface, and B to 
denote the function, so that q = B{fi) on the star. Since g is a constant along streamlines, 
the relationship between the integration variables can be derived from the identity 

dA dzu dB ^^^-j 
dw dfj, dfi 

Next we note that the derivative of the function A{w) takes the form 

where we have used the orthogonality properties of Vg and Vp, and where F = ^ sin 6 (see 
Result 1 from Appendix |A]); note that F = w on the disk surface. Similarly, the derivative 
of the second function B{fi) takes the form 

f - (-sf^ * ■ r (-^ I) . ^ - ■ ^ - ■ • 

where ^ = r/R^ = 1 on the stellar surface. Combining the last three equations allows us to 
write 

U7du7{h;')^ = if.p)^{h;')jfi. (C5) 

This expression allows for a change of variables in either of the two integrals of equation 
(ICll) . and thus shows that the two integrals are equivalent. 
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